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ABSTRACT 



Oh, Context. Many cosmic dust species, among them refractory oxides, form at temperatures higher than 300 K. Nevertheless, most 
J,^ ■ astrophysical studies are based on the room-temperature optical constants of solids, such as corundum and spinel. A more 
» I ' realistic approach is needed for these materials, especially in the context of modeling late-type stars. 

"^ , Aims. We aimed at deriving sets of optical constants of selected, astrophysically relevant oxide dust species with high melting 
j^ ' points. 

Methods. A high-temperature, high-pressure cell and a Fourier-transform spectrometer were used to measure reflectance spectra 
of polished samples. For corundum (Q-AI2O3), spinel (MgAl204), and a-quartz (Si02), temperature-dependent optical constants 
were measured from 300 K up to more than 900 K. Small particle spectra were also calculated from these data. 
Results. All three examined oxides show a significant temperature dependence of their mid-IR bands. For the case of corundum, 
we find that the 13 nm emission feature - seen in the IR spectra of many AGB stars - can very well be assigned to this mineral 
r — ' species. The best fit of the feature is achieved with oblate corundum grains at (mean) temperatures around 550 K. Spinel remains 
^-H , a viable carrier of the 13 /i feature as well, but only for T<300K and nearly spherical grain shapes. Under such circumstances, 
spinel grains may also account for the 31.8 /.tm band that is frequently seen in sources of the 13 /xm feature and which has not 
yet been identified with certainty. 
Conclusions. 
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H 1. Introduction lar composition (e.g. Henning et al. 1999, Henning 2010). 

- - - These databases continue to be indispensable sources of 

Within the past two decades, astromineralogy - the sci- information on the UV to far-infrared properties of min- 

ence of solids in space - has progressed significantly, erals and amorphous solids 

Satellites such as the Infrared Space Observatory (ISO), r^i , a ^ ■^-u r i- -a ■ 

^ J \ /' j^i^g present paper deals with refractory oxides, i.e. 

Spitzer, and Herschel have made it possible to identify . , . , .,, , . 1 u- • + ^^ -k 

-^ ' f J oxide minerals with high melting points. More specih- 

emission and absorption bands of crystalline solids in ,, ,• ^.i, -j • r j + <■ a f 

'^ •' cally, we discuss the mid-mtrared spectra or corundum (a- 

comets, in protoplanetary disks, in the shells of evolved ^^^q^)^ ^^^^^^ (MgAl204), and «-quartz (^-81020 mea- 

stars, in supernova remnants, and in other environments ^^^^^ ^^ ^^^^ temperatures (up to 973 K). Corundum, 

including quasars (cf. Molster, Waters & Kemper 2010). ^^^^^^^ ^^^ q^^^^^ ^^^^ repeatedly been proposed as carri- 

At the same time, the theory of formation and evolution ^^^ ^^ j^^^^^,^^ ^^^^ ^^^^^ ^^^^ .^ astronomical objects (see 
of minerals in accretion disks and stellar outflows has con- 

tributed a lot to our understanding of the gas-solid tran- 1 Note that SiOa and its polymorphs are sometimes classified 

sition m the cosmos (e.g., Gail 2010). Systematic studies ^s silicates because of their structural similarity to tectosili- 

of analog materials in terrestrial laboratories made it pos- cates (see e.g. Klein, Hurlbut & Dana 1999). Nevertheless, we 

sible to establish databases of optical constants of solids will follow the widespread classification of quartz as an oxide 

that are expected or known to form from gases with a so- here, following Strunz & Nickel (2001). 
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below and Sect. 5). Owing to the lack of more comprehen- 
sive data, room temperature dielectric functions, and/or 
powder transmission, spectra of these minerals have been 
used in most previous papers. However, the formation of 
refractory oxides at high temperatures - and the fact that 
these oxides will partly radiate at high temperatures and 
leave their spectral fingerprints in their high temperature 
state - led us to the conclusion that their dielectric func- 
tions need to be derived by in situ- high-T-measurements. 
This term is used to denote measurements of materials in 
their hot state, which is not the same as annealing exper- 
iments. In the latter case, a sample is heated to a certain 
temperature for a while. It may also have undergone an 
irreversible structural transition, and is afterwards cooled 
down to room temperature, at which point the spectra are 
then measured (see, e.g., Jager et al. 2011 for this distinc- 
tion). 

Thus, what we present here - and for the first time 
in a systematic way - are in situ high-T-measurements of 
the IR spectra of a-Al203, MgAl204, and a-Si02, taking 
the effects of anisotropy into account where necessa ry. 
Results of annealing experiments with spinel can be found 
in Fabian et al. (2001). The optical constants derived in 
the present paper will be made publicly available via: 
http : //www . astro . uni- j ena . de/Laboratory/ 
Database/databases . html . 

This paper is structured in the following way: Section 
2 describes the experimental methods and the samples 
that have been used. In Sect. 3, we derive temperature- 
dependent oscillator parameters from our reflectance mea- 
surements. Section 4 presents small-particle spectra for 
different grain shapes. In the final Sect. 5, we discuss as- 
trophysical applications of our data, focusing on the 13 /im 
emission feature in AGB stars. 
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Fig. 1. Sketch of the set up of the HTHP cell in reflection 
mode in the aluminum tank. 



is used (see Fig.[T]). The incident angle of the incoming 
beam on the surface of a sample in the cell is ^ 12°. We 
estimate the relative error in reflectance to the case of per- 
pendicular irradiation to be w 3%. The cell itself has a 
closed sample chamber. The sample holder inside the cell 
is covered by the heater and has a diameter of 13 mm and 
a length of ~15mm. Special rings can be screwed into the 
holder to carry and/or to fix the sample. 

The HTHP cell is supplied by an external temperature 
controller by which the heating temperature and the tem- 
perature gradient can be set. The controller measures the 
temperature of the sample holder and the body of the cell. 
We additionally measured the temperature of the sample 
holder with an NiCr-Ni thermocouple. The temperatures, 
thus externally measured, are only by a few percentage 
points lower than the temperatures measured by the cell. 
All sample temperatures in this paper are externally mea- 
sured ones. 



2. Experimental methods 

2.1. The high-temperature, high-pressure cell 

For all our measurements, we used an IR Fourier- 
transform-spectrometer (Bruker 113v) equipped with a 
water-cooled, high-temperature, high-pressure (HTHP) 
cell. With this cell (Specac P/N 585 0), it is possible to 
heat samples to temperatures up to 1073 K and to take IR 
spectra at the same time. Because of its dimensions, we 
had to deploy the cell into an aluminum tank that is placed 
in one of the two sample chambers of the Bruker spectrom- 
eter. The tank separates the cell environment, which is 
flooded with argon gas at atmospheric pressure, from the 
vacuum inside the spectrometer. The IR beam still can be 
applied for measurements by using two opposite windows 
in the tank. We used both potassium bromide (KBr) and 
polyethylene (PE) windows for measurements in the mid- 
IR (3-25 /xm) and far-IR (25-50 /im) wavelength ranges, 
respectively. 

The HTHP cell, set in reflectance mode, has been 
placed into the tank, where a base plate with two ad- 
justable mirrors for the incoming and outgoing IR beam 



2.2. Samples and sample preparation 

Understanding the changes taking place in refractory 
oxides at high temperatures requires basic information 
on their internal structure. Therefore, we give crystallo- 
graphic information on the respective solids at the begin- 
ning of each of the following sections, starting with AI2O3, 
the most refractory mineral treated here. 



2.2.1. Corundum 

The distribution of Al cations in either a face-centered cu- 
bic or a hexagonal close-packed lattice of oxygen anions re- 
sults in a manifold of polymorphs of AI2O3. When heated, 
all polymorphs undergo phase transitions toward the ther- 
modynamically most stable form, namely corundum (also 
a-Al2 03; Levin & Brandon ll998p . Corundum has a rhom- 
bohedral lattice structure. The oxygen anions are arranged 
in a (slightly distorted) hexagonal close-packing in which 
the aluminum cations occupy two thirds of the octahedral 
interstices. Due to the crystal structure, the interaction 
with light is not isotropic. The optical constants n and 
k (real and imaginary part of the complex refractive in- 
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dex) have to be calculated for polarizations perpendicu- 
lar (ordinary) and parallel (extraordinary) to the crystal's 
symmetry (c) axis. 



Corundum is a highly refractory material. Its con- 
densation temperature is much higher than for most 
other minerals of astrophysical relevance. At pressures of 
~10~^ bar which occur in stellar outflows of evolved stars, 
corundum is already stable at more than 1400 K, while sil- 
icates condense primarily below HOOK (Gail l^OTU)) . As a 
result, corundum should be one of the first condensates to 
form in the envelopes of M type stars. However, quantum 
mechanical calculations of the condensation rates have 
shown that corundum cannot form under these conditions 
by homogeneous condensation from the gas phase. Gail & 
Sedlmayr fl998) point out that the formation of AI2O3 
grains is more likely to occur via grain-surface-reactions 
on already existing grain seeds of other high refractory 
materials such as Ti02 or Zr02. 



For the measurements discussed in this paper, we used 
a plane-parallel disk of pure synthetic corundum, which 
was manufactured by Thorlabs. This disk has polished sur- 
faces and a diameter of 12.7mm, whereas the thickness is 
3 mm. Temperature-dependent reflectivity measurements 
on such pure corundum samples have been presented by 
Gervais & Piriou (|1974ap for temperatures up to 1775 K. 
In their paper, Gervais & Piriou present a description of 
the four-parameter semi quantum (FPSQ) model and a 
comparison with the classical oscillator model. They show 
that the FPSQ model is a more effective way to fit the 
reflectivity of polar crystals in the case of wide reflection 
bands where equal damping of the longitudinal optical 
(LO) and transverse optical (TO) mode cannot be ex- 
pected any more. On the other hand, Thomas et al. (|1998p 
developed a model that does not need a four-parameter de- 
scription of the dielectric function but is based on adding 
temperature-dependent multiphonon contributions to a 
classical oscillator model. Thomas et al. list the param- 
eters for the oscillators, together with their temperature 
dependence, corresponding to a second-order polynomial 
fit. They compared their model with measured data from 
several other authors and the FPSQ model. An applica- 
tion of the results of Thomas et al. to measured emission 
spectra of corundum in the wavelength range of 2.5-20 fim 
is presented in the work of Sova at al. (|1998p . However, the 
presented measurements on corundum in all of these pa- 
pers refer only to the ordinary ray. To the authors' knowl- 
edge, high-temperature IR reflection spectra for the ex- 
traordinary ray that are necessary to calculate small par- 
ticle spectra have not been measured yet. In this paper we 
present for the first time high-temperature IR reflection 
measurements of corundum for the ordinary and extraor- 
dinary ray and the calculate optical constants based on 
the FPSQ model. 



2.2.2. Spinel 

In a wider sense, spinels are minerals with the general com- 
position {Ai^x'Bx)[A.x^2~x]04. The atoms in parenthe- 
sis occupy the tetrahedral sites, and the atoms in square 
brackets occupy the octahedral sites of a cubic close- 
packed oxygen lattice. The sum formula of pure Mg-Al- 
spinel is MgAl204. We will only refer to a material with 
this stoichiometry as spinel in the following. Spinel crys- 
tals with this composition have an elementary cell con- 
sisting of 32 oxygen anions, 16 aluminum and 8 magnesium 
cations. In a perfect spinel crystal, the trivalent Al cations 
are octahedrally coordinated by oxygen ions, while the bi- 
valent Mg cations are tetrahedrally coordinated. However, 
it has been noted by Tropf & Thomas (|199ip and shown 
in detail by Fabian et al. (|200ip that annealing induces 
the Mg and Al ions to change their sites, such that a part 
of the Al ions is occupying the tetrahedral sites and a 
part of the Mg ions is obversely located in the octahe- 
dral sites. This structural transition takes place at 1023- 
1073 K and is apparently irreversible at least at short 
time scales (Tropf & Thomas UMT\ Fabian et al. [2001]) . 
Therefore, synthetic spinel crystals are in particular char- 
acterized by this deviation from the perfect crystal, since 
high temperature conditions are needed for their produc- 
tion. In contrast, natural spinel crystals that were exposed 
to low-temperature conditions for long periods show an 
almost perfect spinel crystal structure. Apparently, under 
low temperatures and after long periods, the aluminum 
cations leave the occupied tetrahedral sites and move to 
the octahedral sites, while the opposite is true for the mag- 
nesium cations (Fabian et al. I2001|) . However, in the lim- 
iting case with zero Mg content, three eighths of the Al 
cations occupy tetrahedral sites, and the spinel transforms 
into 7-AI2O3. 

Spinel has a very high melting point of ~2300K at at- 
mospheric pressure (Tropf & Thomas ll99ip - even higher 
than for terrestrial corundum. However, at the pressures 
prevailing in the dust-formation zones of AGB-stars (in 
the range of ^10^^ bar), corundum is able to form first 
condensates at higher temperatures (at ~1400K) than 
spinel (which we expect to condense below 1150 K). As 
mentioned above, spinel has a cubic crystal structure 
and is therefore optically isotropic. In our measurements 
we used a (pure) synthetic spinel disk (manufactured by 
SurfaceNet) with a diameter of 12.7mm and a thickness of 
0.1 mm. Like the corundum disk the spinel sample has pol- 
ished surfaces. High-temperature data have already been 
obtained by Thomas et al. (|1988p and by Sova et al. 
(|1998p . Thomas et al. published absorption coefficients 
in the range of 2.5-10 /xm for various temperatures up to 
2000 K and a comparison of the results with the multi- 
phonon model. Sova et al. measured the emissivity of a 
pure spinel disk over a wavelength range of 2-12.5 /im at 
~2000K and compared their results with the combined 
model of multiphonon contributions to the classical dielec- 
tric function from Thomas et al. (|1998p . Our data cover 
the whole wavelength range from 5 /zm to 50 /im. 
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2.2.3. Quartz 

Quartz is one of the three low-pressure phases of crys- 
talline Si02- The other two are tridymite and cristobaliteo 
It can be divided into a lower temperature, lower symme- 
try (a) and a higher temperature, higher symmetry (/3) 
modification. 

The mineral a-quartz has a trigonal lattice structure 
at which the unit cell includes three silicon and six oxy- 
gen ions. The silicon cations are tetrahedrally coordi- 
nated by four oxygen anions which creates the [Si04]^~- 
tetrahedron, the basic structure of almost all silicates 
and silica polymorphs. Each oxygen anion belongs to two 
tetrahedra, which leads to a three-dimensional grid of 
interconnected tetrahedra, with triples of them winding 
about a hypothetical axis that is parallel to the c-axis. 
Thereby, a helix-like structure is formed. It should be 
mentioned that these structures have no particular influ- 
ence on the chemical and physical behavior of the crystal, 
they are just a help in visualizing the basic structure of 
a-quartz. 

The change between a- and /3-quartz is taking place at 
846-847 K at atmospheric pressure by a displacive phase 
transition (Lakshtanov et al l2007p . /3-quartz has a hexago- 
nal lattice structure and its tetrahedra are arrangeded al- 
most in the same way as for a-quartz. The only difference 
lies in the alignment of the tetrahedra. In this respect, /3- 
quartz has a higher symmetry than a-quartz. By cooling 
down a /3-quartz crystal, the opposite transition occurs, 
and an a-quartz is formed again. 

Quartz has been the topic of many studies, since it is 
one of the most important minerals on Earth and, with 
its physical and chemical properties, an outstanding ma- 
terial in science and technological applications. Spitzer 
& Kleinman (|196ip were the first to examine the IR 
reflection properties of a-quartz with respect to classi- 
cal dispersion theory. They presented very comprehen- 
sive room-temperature studies of natural and synthetic 
quartz crystals. Gervais & Piriou (|1975|) presented the 
first temperature-dependent IR reflection data of quartz. 
They fit their data between 295 and 975 K with the al- 
ready mentioned FPSQ model (Gervais & Piriou I1974ap 
and analyzed the changes in the band parameters with 
temperature, especially in the range of the phase transi- 
tion from a~ to /3-quartz. 

In the present study, we took IR reflection measure- 
ments on a natural a-quartz crystal from Brazil (see Fig. 
[5]). The crystal has a maximum length, height, and width 
of II, 7.5, and 8mm, respectively. EDX analysis did not 
reveal any case of impurities or deviations from stoichiom- 
etry. Because of the crystal structure of a-quartz, the 
interaction with light is anisotropic. As in the case of 
corundum, one has to distinguish between polarizations 
parallel and perpendicular to the c-axis. Fortunately, be- 
cause of the distinct hexagonal shape of the natural crystal 




The high-pressure phases of Si02 are coesite, stishovite, 
and seifertite. In total there are ten Si02 polymorphs (Klein, 
Hurlbut & Dana 1999; El Goresy at al. 2008). 



Fig. 2. Photograph of the quartz crystal that has been 
used for the measurements. The diameter of the field of 
view is 25 mm. 



used, we could locate the direction of the crystallographic 
c-axis quite easily. Moreover, the crystallographic c-axis 
lies parallel to the large natural grown surfaces of our a- 
quartz crystal. Thus, spectroscopic measurements parallel 
and perpendicular to the c-axis could be taken on one of 
these surfaces. The crystal was cut (10 mm x 5 mm) so 
that it fit into the HTHP cell where the surface used for 
the measurements has been polished with diamond paste 
of 0.25 /^m fineness. 

2.3. Reflectance measurements 

We performed IR measurements with the HTHP cell in 
reflection mode (see Sect. 12. ip in a wavelength range of 
5-50 /im. As mentioned above, special rings were used to 
hold the samples. There are two kinds of these rings: one 
with a hole of 12 mm that was used to hold the synthetic 
samples and one with a hole of 5 mm for the quartz crys- 
tal. For corundum and quartz, a polarizer was placed in 
the outgoing beam to distinguish between the different 
crystal orientations. We used two kinds of polarizers: a 
KRS-5 polarizer for mid-IR measurements (5-25 /im) and 
a PE polarizer for far-IR measurements (25-50 /im). With 
the aid of a low-resolution, real-time spectra shooting in 
the setup mode of the spectrometer software, the correct 
polarizer adjustment was established by finding the max- 
imum of bands or by the disappearance of other bands. 

As reference for each measurement, a spectrum of a 
gold mirror with a diameter of 12.7 mm, taken at room 
temperature (300 K), was used. Additionally, IR reflec- 
tion spectra were taken after the measurements from the 
empty cell only equipped with the ring that was used for 
the respective sample. These empty measurements were 
also taken at the respective temperature of the measured 
sample and were necessary to exclude influences of the 
whole sample holder setup. The spectra of sample Rs and 
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empty cell Re, both related to the reference, define the 
pure reflection spectra (i?) via 



R = 



Rs — Re 
I -Re 



(1) 



It could be verified that the heated samples and the 
sample holders themselves serve as additional IR sources 
whose emissions contribute to the spectra. While the ra- 
diation directly emitted toward the detectors gives a DC 
signal and is filtered out, the radiation back to the inter- 
ferometer of the spectrometer becomes modulated and is 
superimposed on the modulated IR beam from the globar 
source of the spectrometer. Therefore, all interferograms 
would actually be superpositions of two sources. To get 
rid of the additional influence, we measured the interfer- 
ogram of the heated cell alone by decoupling the globar 
source and subtracted the result from each interferogram 
with activated globar source. After this, the resulting in- 
terferogram is Fourier-transformed into the desired (single 
channel) reflection spectrum that is used to calculate Rs 
or Re- 

For the heating process, we used a constant temper- 
ature gradient of 10 K/min to reach the respective mea- 
surement temperature. For corundum and spinel, we took 
mid-IR spectra at 300 K, 551 K, 738 K, and 928 K. For 
quartz, we also measured at 833 K to derive more infor- 
mation on the development of the spectrum in the vicinity 
of the phase transition from a— to /3-quartz. We chose no 
more than these individual temperatures to prevent the 
heater from stronger abrasion. 

3. Derived data and data analysis 

3.1. Derivation of the optical constants 

To calculate the optical constants from reflection measure- 
ments, we fitted our data with the FPSQ model of the 
dielectric function. This model provides the opportunity 
to fit measured data with independent damping param- 
eters for the LO and the TO mode, while the classical 
oscillator model only allows to assign different frequencies 
to the above mentioned two modes. The necessity of as- 
suming independent damping parameters arises in crystals 
with a strong mode splitting, where LO and TO phonon 
modes should have different phonon decay times (Gervais 
& FirioulMig). 

Based on the classical dispersion theory, the dielectric 
function can be written as: 






A< 



r2? 



^n] 



-OJ^ 



IJjL 



(2) 



where Ae^, flj, and jj are the strength, resonance fre- 
quency, and damping of the jth oscillator. According to 
Maxwell's equations, the transverse modes are the poles 
of the dielectric function, and in contrast the longitudi- 
nal modes are the zeros. Barker (|1964p has shown that by 



multiplying all resonant denominators to Equation ([2]), a 
polynomial of degree 2n is created (n is the number of os- 
cillators) and that this polynomial can be factorized with 
the zeros of the dielectric function. Consequently, the di- 
electric function is extended with the LO frequency and 
damping parameters, and it can be rewritten in its alter- 
native form: 



ri2 



^ - njLo^ 






(3) 



This is the FPSQ model in its simplest representation. 
A deeper treatment of this problem leads to the quantum 
form of the dielectric function (Gervais & Piriou I1974a|) 
that goes beyond the classical model toward a quantum 
mechanical description of the dispersion theory of anhar- 
monic crystals. In our relatively simple case of fitting re- 
flection data at high temperatures, the quantum form will 
not be important, though it becomes important in dis- 
cussing the temperature dependence of the mode frequen- 
cies and the dampings, so it is rudimentarily discussed in 
Sect.O 

The optical constants n and k (the real and the imag- 
inary part of the refraction index) are a function of e rep- 
resented by the formulae: 



/? + P 



2n 



cl'2 



(4) 



(5) 



The reflection R on a surface at near normal incidence in 
a vacuum can be expressed as a function of n and k: 



R 



(1 



(1 



p- 



(6) 



Together with Eqs. (O, Q, and ([S]) measured IR reflection 
spectra can be fit with the five parameters of this formula. 
For all data, the fit procedure has been performed by 
using the datafit-innction of the software SCILAB 5.4.0 
(SCILAB consortium [Mil)- 

3.1.1. Crystalline a-Al203 (corundum) 

In Figs. [3] and |4] the measured reflectivities of the corun- 
dum sample for the ordinary and the extraordinary ray at 
four different temperatures are shown, together with their 
best fit, with the aid of the FPSQ- model (Eq. ©). As 
can be seen, the main changes due to heating occur in the 
band positions and the intensity of the bands. In general, 
the bands tend to shift toward longer wavelengths, while 
the intensity decreases with increasing temperature. This 
effect is due to the anharmonicity of the atomic potentials 
and can be seen more clearly in the values of the oscillator 
parameters listed in Tables [1] and [51 Here, the mode po- 
sitions move toward longer wavelengths (lower wavenum- 
bers, respectively), and the respective damping constants 
increase with increasing temperature. 
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Fig. 3. Comparison of the fittings and the reflection spec- 
tra for the ordinary ray (i^u-type modes) of the corundum 
sample at four different temperatures. The fittings have 
been done with four oscillators. 



Fig. 4. Comparison of the fittings and the refiection spec- 
tra for the extraordinary ray (v42ti-type modes) of the 
corundum sample at four different temperatures. The fit- 
tings have been done with three oscillators. 



According to the prediction of the irreducible repre- 
sentation of the vibration modes of corundum flishi ll975)) . 
the spectra for the ordinary ray have been fit with four os- 
cillators (£'„-type modes). In contrast, spectra for the ex- 
traordinary ray have been fitted with three oscillators, al- 
though theory predicts only two oscillators to be IR active 
(^2ji-type modes). Though, the third oscillator is needed 
to fit the distinct sharp minimum in the reflection spectra 
around 21 /im. An additional band in that range might 
indicate effects of surface polishing or of imperfect polar- 
ization. 

The quality of the fits is different, depending on tem- 
perature and wavelength range. In general, the noise and 
the infiuence of water vapor in the FIR is much stronger 
than in the MIR, hence the fitting performance, especially 
for weak bands, is lower, and the quality of the result- 
ing fit decreases for higher temperatures. As can be also 
seen in Figs. [3] and |4l the maxima of the measured re- 
flectance spectra were not flt exactly. The measured val- 
ues are either somewhat too high (for _E„-type modes) or 
too low (for A2„-type modes). Barker (|1964p has shown 
that the kind of surface treatment can have a strong in- 
fluence on the measured reflectance spectra proflles owing 
to additionally excited forbidden modes. Especially the 
structures on top of the strongest bands measured in both 
polarization directions give rise to the statement that the 
differences between measurement and fit in that range are 
caused by these forbidden modes. Likewise, these struc- 
tures become weaker with increasing temperature. 

3.1.2. Crystalline MgAl204 (spinel) 

In Fig. \5\ the spinel refiection spectra (dotted lines) are 
compared with their best fit (solid lines) at four different 
temperatures. The fits are in very good agreement with 



Table 1. Comparison of the fitting parameters of corun- 
dum for the Eu-type modes (E_Lc) at different tempera- 
tures. The definition of the parameters refers to Eq.([3l). 



j r(K) Utoj 
{cm-'-) 



iTOj 



{cm~^) {cm~^) 



1 300 


384.02 


6.03 


387.46 


5.18 


3.05 


551 


380.83 


6.13 


383.22 


6.55 


3.06 


738 


379.9 


9.86 


382.36 


8.56 


3.08 


928 


377.51 


10.97 


379.8 


10.1 


3.1 


2 300 


439.22 


3.23 


481.96 


2.96 


_ 


551 


436 


7.51 


477.49 


3.86 




738 


432.94 


8.19 


474.39 


6.99 




928 


430.21 


10.49 


471.18 


10.15 




3 300 


569.35 


7.86 


908.23 


22.37 


_ 


551 


564.42 


8.89 


901.64 


27.84 




738 


560.1 


10.24 


898.5 


36.78 




928 


555.46 


14.82 


894.5 


45.36 




4 300 


634.36 


5.6 


630.59 


8.53 


_ 


551 


628.8 


8.13 


624.95 


11.37 




738 


624.51 


10.8 


620.9 


14.32 




928 


619.99 


14.21 


616.66 


17.55 





the measured data and can be treated as almost perfect 
models for the measurements within the accuracy given 
by the line thickness. 

For a normal, non-defective MgAl204 crystal, group 
theory predicts only four IR active oscillators whose TO- 
modes are located at around 15, 18, 21, and 33/Ltm 
f Thibaudeau 12002]) . However, the weak band at 43 /xm, as 
well as the weak shoulders at around 12 and 19 /im, made 
it necessary to use eight oscillators for the fitting proce- 
dure (see Table ^ . These additional modes are forbidden 
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Table 2. Comparison of the fitting parameters of corun- 
dmn for the A2u-iyp& modes (E||c) at different tempera- 
tures {j = 1,2). The third oscihator could not be verified 
as an A2«-type mode. 



1000 



wavenumber {cm'') 
500 400 



300 



200 



by group theory but can be excited by the mentioned dis- 
ordered cation distribution in synthetic spinels (Fabian et 
al. (|200ip ') and by defects in the crystal lattice similar to 
the case of corundum (see above). 

Previous studies have also made use of additional 
modes by fitting reflection spectra of spinel. Fabian et al. 
(|200ip used a classical Lorentz oscillator model and de- 
rived fits with eight oscillators in a wavelength range of 
9-35 jum, which means they did not treat the 43 /im band. 
Therefore, they used an extra mode at 11.4 /im (5-mode) 
to fit a weak shoulder on the short-wavelength edge of 
the broad band between 11-15 /^m, which they link to the 
Al/Mg-ratio of their spinel samples and the partial fill- 
ing of tetrahedral interstices by Al- instead of Mg-ions. In 
our spectra we could not see any trace of such a shoul- 
der. Thibaudeau et al. p002)) applied the FPSQ-model 
to their spectra with six oscillators. Based on their re- 
sults, we added two more oscillators to the short- and the 
long- wavelength edges of the broad 19 /im band to fit the 
shoulder at 17 /im and the decrement toward longer wave- 
lengths correctly. These extra oscillators have already been 
found in previous studies (Fabian et al. (|200ip . O'Horo et 
al. ((inZSl), Preudhomme & Tarte PWT]) ) and might be 
due to the ratio of Mg/Al-cation disorder and defects, re- 
spectively. 

Comparing the results of the fitting parameters listed 
in Tab.[3l it is noticeable that only the strongest oscillators 
show a monotone development of their parameters with 
the temperature. The weaker oscillators (j = 1,3,6) do 
not exhibit such distinct behavior, which might be because 
in general the fitting of weak modes strongly depends 
on the quality of the measured spectra, making the fits 
non- unique (e.g. the weak band around 43/zin). However, 
in general the mode positions move toward longer wave- 



1.5 - 



1.25 



j 


T(K) 


ilroj 


7TOj 


^LOj 


ILO] 


too 


1 






(cm-^) 


(cm-^) 


(cm-i) 


(cm-^) 






1 


300 


399.68 


4.68 


511.05 


1.42 


2.9 


^0.75 




551 


394.76 


11.2 


508.09 


2.79 


2.91 






738 


391.91 


18.58 


505.24 


2.88 


2.91 


0.5 




928 


390.54 


22.57 


502.65 


3.22 


2.97 


0.25 


2 


300 


481.58 


3.42 


480.93 


3.21 


- 






551 


476.38 


7.83 


475.7 


6.99 









738 


472 


13.3 


471 


12.13 








928 


471.19 


17.11 


469.84 


16.32 






3 


300 


582 


4.17 


884.75 


21.57 


_ 


Fig. 5 




551 


577.36 


8.93 


882.32 


35.96 






738 


573.5 


11.37 


878.47 


44.71 








928 


573.35 


15.74 


874.39 


49.24 








15 20 25 30 35 40 

wavelength (mn) 



Reflection spectra of the spinel crystal at four dif- 



lengths upon heating, while their respective damping in- 
creases with temperature. 

3.1.3. Crystalline a-Si02 (quartz) 

In contrast to spinel and corundum, the spectra of quartz 
show a lot of distinct narrow bands, as can be seen in 
Figs. inland [7] for the ordinary and the extraordinary rays, 
respectively. Theoretical calculations predict eight E- and 
four A2-type modes to be IR active in a-quartz, corre- 
sponding to the excitations in the ordinary and extraordi- 
nary rays, respectively. By contrast, in /3-quartz, four E- 
and only two A2-type modes should be IR active as ana- 
lyzed by Scott and Porto (|1967p . Two A2 modes remain 
from the a-/3-transition of quartz, while the other two be- 
come forbidden. The same holds true for the E modes, but 
here the missing modes only become IR inactive, while 
they still remain Raman active (Scott and Porto [TM7| . 

The fitting parameters for the ordinary and the ex- 
traordinary ray at five different temperatures are shown 
in Tables S] and [S] We found seven E- and four A2- 
type modes for temperatures below 847 K. The eighth E- 
type mode should be found at wavelengths around 78 /zm, 
which is, however, beyond our covered wavelength range 
and therefore has not been measured. Anyway, the pa- 
rameter values are quite consistent with data given by the 
literature (Scott and Porto [T^7i Gervais & Piriou [TOTSIl 
and they fit our measured data very closely (see Figs. [S] 
andlZl). 

By reaching 928 K for the ordinary ray, it can be clearly 
seen that the bands at 8.2, 14, 25, and 37 /im (mode num- 
ber 7, 4, 2, and 1) disappeared. According to Scott and 
Porto, these are the four modes that should become IR 
inactive by the transition to /3-quartz. Already at 833 K, 
the modes 7, 4, and 1 became very weak so that there 
is almost no evidence of their existence in the spectrum. 
In the extraordinary ray, the spectra also show the tran- 
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Table 3. Comparison of the fitting parameters of spinel. 
Oscillators 2,4,5, and 7 represent the predicted Tiu-type 
modes of the theory. 



wavenumber {cm' ) 
500 400 



T(K) a 



■TOj 



iTOj 



^ 



LOj 






1 


300 


231.05 


8.97 


231.58 


8.29 2.77 




551 


226.96 


9.99 


227.54 


9.41 2.8 




738 


230.02 


9.79 


230.34 


9.18 2.82 




928 


227.33 


13.81 


227.78 


13.45 2.77 


2 


300 


307.12 


12.65 


310.73 


13.68 




551 


304.4 


16.38 


308.36 


17.02 




738 


302.2 


20.15 


306.49 


20.68 




928 


300.07 


20.15 


304.01 


20.58 


3 


300 


427.11 


134.26 


442.85 


112.26 




551 


426.34 


139.54 


438.44 


108.65 




738 


438.12 


140.41 


444.8 


96.33 




928 


440.49 


140.09 


442.56 


97.23 


4 


300 


481.53 


30.58 


608.03 


37.88 




551 


475.24 


35.89 


598.4 


42.79 




738 


469.66 


38.7 


590.55 


48.62 




928 


466.52 


45.67 


582.36 


55.59 


5 


300 


531.37 


55.1 


558.16 


70.14 




551 


527.69 


63.05 


553.85 


74.68 




738 


521.27 


65.86 


549.66 


79.03 




928 


520.32 


76.53 


546.93 


82.07 


6 


300 


581.8 


74.64 


546.32 


120.96 




551 


576.95 


74.12 


547.55 


122.79 




738 


569.98 


75.11 


542.25 


126.64 




928 


563.9 


73.29 


544.39 


126 


7 


300 


672.28 


36.67 


866.53 


39.52 




551 


662.28 


47.55 


862.97 


47.55 




738 


655.42 


55.98 


860.15 


55.98 




928 


649.32 


64.24 


855.98 


64.24 


8 


300 


811.63 


77.37 


800.67 


73.77 




551 


797.84 


86.36 


785.72 


76.12 




738 


789.47 


92.02 


776.78 


79.87 




928 


786.4 


93.61 


772.56 


85.91 



sition to /3-quartz. Here, modes number 1 and 3 (at 26 
and 13 /xm) disappeared by reaching 928 K, as predicted 
by theory. These results correspond very well to the data 
taken by Gervais & Piriou (|1975|) . who comprehensively 
investigated the phase transition of quartz with IR spec- 
troscopy between 7.7 and 33.3 /im. At this point we cab 
complement their study with data on the temperature de- 
velopment of the E-type mode at 37 jum that has not been 
treated in their work. 

The quality of the fits decreases at longer wavelengths 
due to noise and absorption bands of water vapor in 
the measured spectra. In those with very strong water- 
absorption bands, we tried to keep the fits slightly above 
the measured spectra at longer wavelengths since the sig- 
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1.25 - 




0.75 
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Fig. 6. Reflection spectra of the quartz crystal in the or- 
dinary ray at five different temperatures {E-type modes). 



wavenumber {cm'') 
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Fig. 7. Reflection spectra of the quartz crystal in the ex- 
traordinary ray at five different temperatures (yl2-type 
modes) . 



nal intensity had decreased in that range (see e.g. Fig. 

El. 
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Table 4. Comparison of the fitting parameters of quartz 
for the ordinary ray. 



Table 5. Comparison of the fitting parameters of quartz 
for the extraordinary ray. 



j T(K) n 



TOj 



iTOj 



^ 



LOj 



iLOj 



{cm ^) {cm ^) {cm ^) {cm ^) 



1 300 


264.13 


6.68 


265.36 


7.08 


2.34 


551 


257.69 


9.21 


258.36 


9.05 


2.36 


738 


255.06 


9.29 


255.77 


9.87 


2.38 


833 


256.4 


11.18 


256.69 


10.94 


2.39 


928 


- 


- 


- 


- 


2.36 


2 300 


393.55 


4.25 


402.88 


4.33 


_ 


551 


395.03 


5.26 


401.96 


4.88 




738 


396.5 


10.44 


402.99 


11.05 




833 


400.14 


12.7 


404.6 


11.83 




3 300 


449.78 


3.43 


507.89 


2.87 


_ 


551 


445.61 


4.9 


504.58 


7.88 




738 


438.27 


9.28 


499.84 


9.29 




833 


432.7 


11.19 


495.08 


10.62 




928 


423.04 


11.85 


491.45 


12.18 




4 300 


694.3 


11.91 


697.31 


12.11 


_ 


551 


695.37 


13.54 


697.13 


12.52 




738 


688.59 


15.17 


688.99 


14.51 




833 


680.64 


38.36 


680.51 


37.27 




5 300 


795.36 


7.75 


808.72 


8.11 


_ 


551 


790.2 


13.44 


804.58 


13.02 




738 


786.14 


16.97 


799.52 


16.62 




833 


782.98 


19.92 


796.81 


20.2 




928 


780.28 


22.85 


793.68 


26.32 




6 300 


1064.57 


6.98 


1226.34 


10.91 


_ 


551 


1063.3 


9.31 


1221.61 


18.66 




738 


1058.36 


15.17 


1217.97 


20.43 




833 


1057.85 


20.31 


1220.69 


31.12 




928 


1060.37 


19.53 


1226.38 


43.43 




7 300 


1160.03 


8.38 


1157.02 


7.98 


_ 


551 


1157.62 


10.57 


1155.85 


9.98 




738 


1155.9 


36.83 


1153.29 


34.45 




833 


1154.85 


41.18 


1154.19 


39.93 





r(K) ^TOj 

{cm-^) 



{cm-^) 



n 



LOj 



{cm 



JLO] 

{cm-^) 



1 300 


363.03 


5.18 


386.42 


3.76 


2.39 


551 


369.72 


7.89 


388.87 


11.92 


2.38 


738 


376.11 


15.02 


392.48 


18.9 


2.38 


833 


386.59 


20.1 


396.64 


20.05 


2.38 


928 


- 


- 


- 


- 


2.44 


2 300 


494.99 


10 


552.82 


5.3 


_ 


551 


486.1 


17.85 


547.22 


10.28 




738 


473.93 


23.83 


540.15 


14.12 




833 


465.2 


24.9 


533.84 


17.53 




928 


442.2 


30.04 


527.43 


24.19 




3 300 


775.58 


8.48 


789.17 


6.51 


_ 


551 


775.09 


13.78 


785.37 


15.29 




738 


775.27 


19.98 


782.47 


20.58 




833 


777.02 


31.45 


782.29 


30.85 




4 300 


1073.37 


15.21 


1238.42 


16.41 


_ 


551 


1070.08 


20.03 


1237.02 


30.04 




738 


1066.91 


22.76 


1234.71 


39.73 




833 


1065.91 


26.52 


1239.85 


50.93 




928 


1064.14 


32.9 


1240.89 


55.58 





ujj and a frequency shift Aujpij representing the influ- 
ence of the phonon interactions. According to Gervais and 
Piriou (|1974al b). flj takes the form 



m 



2ll!jAujpij. 



(7) 



Owing to the anharmonicity, both Uj and Aajp/j|f| 
hence ftj, show a temperature dependence. While for 
Aujpjj the temperature dependence arose from the mode- 
specific excitation probability of the phonons {Bose- 
iJmstem-distribution) , the temperature dependence of ujj 
is due to the thermal expansion of the solid and can be 
approximated by the model of Grtineisen (Lowndes 119701 
Jasperse et al. 119661 Gervais and Piriou I1974ap with: 



3.2. The temperature dependence of the oscillator 
parameters 

3.2.1. Theoretical point of view 

The lattice potential energy of a solid is, in general, an- 
harmonic (Maradudin and Fein 119621 Cowley [1963]) . The 
anharmonicity has, in principle, two main effects on the 
solid and its spectral properties: thermal expansion of the 
solid and interactions of phonons of different modes in- 
fluence their energy and give them a finite lifetime (the 
reciprocal counterpart to the damping) . If the phonon in- 
teractions are sufficiently small, the phonon frequency Qj 
for the jth mode as introduced in ^ can be expressed in 
terms of an unperturbed quasiharmonic mode frequency 



Wj = iooj + Sujj = ujQj — loqj / gj{T')a{T')dT\ (8) 

Jo 

where ujqj represents the phonon-mode frequency at 
OK for the jth mode. Here, gj is the respective mode 
Griineisen parameter and a the volume thermal expan- 
sion of the material. While for most of the materials a 
can be approximated well by a linear temperature depen- 
dence, gj has a rather complex relation to the temperature 
(Jasperse et al. I1966p . which is often unknown. However, 



^ It should be noted that Aujpij in general carries a fre- 
quency dependence that still can be neglected in the FPSQ 
model (Gervais and Piriou [T974ap . Also, ilj is written here as 
representative of both LO and TO mode frequencies. 
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in many cases gj can be assumed as a constant value rang- 
ing between and 2 (Lowndes 119701 Gervais and Piriou 
I1975[l . in which case 6ujj will show a quadratic tempera- 
ture dependence. 

The phonon-interaction contribution to the frequency 
shift Aujpij has been determined by Maradudin and Fein 
(|1962|) and by Cowley (jl963p involving three phonon pro- 
cesses to second-order and four phonon processes to first- 
order pertubation theory. With their assumptions, they 
obtain a linear temperature dependence of the frequency 
shift in the high-temperature limit. A quadratic tempera- 
ture dependence can be obtained by including terms of 
higher order pertubation theory and/or higher phonon 
processes as shown by Ipatova et al. (|1967p and Gervais 
and Piriou (jl974bp . By assuming Aujpij to be sufficiently 
small, equation ([7]) can be simplified to 



flj — UlQj + Su!j + AuJpjj . 



(9) 



Maradudin and Fein (|1962l) pointed out that the fre- 
quency shift Aujpij appears as the real part of a complex 
phonon self energy. In their theory, the imaginary part 
is represented by a function Tj that has the property of 
altering the phonon energy, so it can be seen as a repre- 
sentative of the damping function. Like Aujpij , damping 
appears as a consequence of the phonon interaction due 
to anharmonicity, specifically as the inverse of a now finite 
phonon lifetime. The quantity Tj is related to the classical 
damping factor jj by the equation 



7j = Tj 



2wi 



'oj 



(10) 



j^Qj 



Like the frequency shift, the damping generally shows 
a quadratic temperature dependence in the high tempera- 
ture limit (Ipatova et al. 119671 Gervais and Piriou [T974bp . 

3.2.2. The case of corundum 

The temperature dependence of the complete set of mode 
parameters for corundum is shown in Fig.[8l For this com- 
parison we chose to plot the relative changes of the pa- 
rameters with respect to their respective values at 300 K, 
similar to the low temperature forsterite data presented 
in Koike et al. (|2006p . A first look shows that the thermal 
shift of the resonance frequencies flro (upper left panel) is 
relatively similar for all modes; i.e., the Uto values change 
by 1.5%-2.5% for 900 K relative to 300 K, while for all 
other parameters describing damping and strength of the 
resonances the differences among the modes are consider- 
ably larger. For some reason, the damping parameters of 
the A2u-type modes increase more than for the Eu-type 
modes, however, no other clear trends (e.g. strong vs. weak 
modes) can be found. 

A fitting of the temperature dependence with respect 
to the given theory is quite difficult because of the limited 
temperature range and resolution, especially without any 
data for lower temperatures. Nonetheless, we try fitting 
both Q.J (T) and jj (T) by the quadratic equations 



^(T) = f7j(300)+^ji(T-300)+/3j2(r2-3002) (11) 

and 

lAT) = 7^(300) + -/jiiT - 300) + 7,2 (T^ - 300^), (12) 

respectively. In these, we have subtracted the values at 
300 K on both sides, leaving the coefficients /3ji, f3j2, 7ji, 
and 7j2 unchanged with respect to a K-based fitting 
curve. The coefficients have been 

restricted such that flj (T) is ensured to monotonically 
decrease, at least until the melting point of corundum 
(ss2300K), while jj{T) has to monotonically increase. 
Thus, Pji should be normally negative, and /3j2 is not 
allowed to be larger than 1 x 10"^ cm~^K~^, whereas 7^1 
and 7j2 are both restricted to positive values. 

The results are listed in the Tables [6] and [T] As can 
be seen, only the last A2«-type mode has a positive /3ji, 
which is however compensated for by a large negative 
/3j2, indicating that the shift of the mode becomes much 
stronger at high temperature. For some of the E^ and most 
of the A2u-type modes, /3j2 reaches the positive limit of 
1 X 10~'^ cm~^K~^, which indicates that the data imply 
a significant weakening of the temperature trend at high 
T. Generally, the absolute value of /3ji tends to increase 
with increasing mode frequency, which is because we fit 
the absolute frequencies here, having already noted above 
(for the TO frequencies) that the behavior of the relative 
frequency changes would be similar among the modes. 

The behavior of the 7j(r) with temperature appears 
to be quite diverse, as indicated by the fitting parameters. 
For the A2u-type modes, where the damping is generally 
rising more sharply with the temperature than for the -E„- 
type modes, the changes can more easily be represented 
by a linear fitting (large 7ji) and the trend tends to be 
flattened at higher T (that would result in a negative 7^2 
but is "zero" due to our restriction). In contrast, for the 
modes in the ordinary ray (Eu-type modes), considerable 
positive quadratric behavior (7^2 dominating) is found. 

Trying to distinguish between thermal expansion ef- 
fects and the phonon-interaction contribution in the fre- 
quency shift, we write Atupij = ApjjiT + Apij2T^ and 
set gj as temperature-independent (Eq.[5]). Then, /3ji and 
/3j2 can be expressed by 



^2 



Apiji — ujQjgjao 
Apij2 - ujojgjai/2, 



(13) 



where ag and ai are the coefficients of the linear ther- 
mal expansion a — ao + aiT. For corundum, they have 
been determined by Piquet et al. (|1999p to ag — 2.081 x 
10~^ and ai = 6.6 x 10~^K~^. At room temperature, 
Lodziana and Parlinski (|2003p give an interval for all mode 
Griineisen parameters of corundum ranging from 0.4 to 
1.8 with an averaged value {g) = 1.14. They find that the 
highest values (1.8) are reached by the A2U and the lowest 
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Fig. 8. Comparison of the temperature development of the mode parameters of corundum as given in Tables [T] and [H 
The parameters have been normalized to their respective value at 300 K and are ordered with respect to TO and LO 
mode, as well as to frequency and damping. 



Table 6. 

modes). 



Comparison of the temperature-fit parameters of corundum for the modes in the ordinary ray (E'u-type 



J 


Al 




Iij2 




7ji 




l32 




(IQ-^cm" 


-1^-1) 


(10-^cm" 


-1^-2) 


(10-^cm- 


-^K-^) 


{IQ-'^cm-^K-^) 


1 {TO) 


-9.84 




0.2 









6.6 


2 {TO) 


-10.67 




-3.1 




12.02 







3 {TO) 


-14.2 




-6.5 









7.9 


4 {TO) 


-19.91 




-2.4 




2.61 




9 


\{LO) 


-13.62 




1 




3.38 




3.7 


2{L0) 


-18.39 




1 









8.9 


3 {LO) 


-23.54 




1 









30 


A{LO) 


-22.72 




0.4 




6.05 




6.8 



(0.4) by the i?„-type modes. We adopt these values for all 
temperatures. 

Taking the room-temperature resonance frequencies as 
an estimate of the wqjj we can derive, for the latter {E^- 
type modes, gj — [0.4,1.14]), a significant Ap/ji in the 



range between —14.63 x 10 "^ and —0.7 x 10 



'K- 



Similar to the mentioned development of /3ji, the abso- 
lute value of Ap/ji also tends to increase with increasing 
rij(300). On the other hand, Apjj2 can be positive or neg- 
ative depending on the mode, and it lies between — 5 x 10~^ 
and 3.09 x 10^^ cm^^i^^^. It shows no clear development 
with increasing r2j(300). 

The estimated values of Ap/^i for the j42u-type modes 
in the limits gj — [1.14, 1.8] are only for j = 1 in 



the same magnitude as the ones for the i?u-type modes. 
In the case of j — 3, they vary between —3.2 x 10~^ 
and 4.8 x 10^^ cm~'^K~^. Nevertheless, Apjj2 is posi- 
tive for both modes, and it lies between 4 x 10^^ and 
7.91 x 10~^ cm~'^K^^. We found an overall increase in the 
values of both Ap/ji and Apij2 with increasing r2j(300). 
Since both mechanisms (phonon interaction and thermal 
expansion) generally have negative contributions in their 
linear coefficients, the calculations confirm that the values 
of /3ji are reasonable. 
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Table 7. Comparison of the temperature fit parameter of corundum for the modes in the extraordinary ray {j — 1,3: 
^2M-type modes). 

j Pjl Iij2 7jl 1]2 

(10-^cm-ii<-i) (lO-^cm-^A--^) (IQ-^cm-^T^-') {IQ-'' cm-^ K'^) 



3.4 



16.3 



1 {TO) 


-17.12 


1 


2 {TO) 


-19.65 


1 


3 {TO) 


-17.02 


1 


l{LO) 


-9.34 


-I 


2{LO) 


-20.59 


1 


3 {LO) 


3.28 


-] 



29.21 





15.48 


5.4 


15.53 


2.1 


3.25 





8.25 


10.5 


47.94 






3.2.3. Spinel and quartz 

For spinel, the temperature fit is much more difficult due 
to the weak bands and the unknown degree of ion disorder 
that can have a big influence on the fits in the case of a 
synthetic spinel (Thibaudeau et al. I2006p . Therefore, be- 
cause of the missing agreement with the mode parameters, 
we will not show the temperature fits of spinel. 

The temperature dependence of the phonon frequen- 
cies and dampings in quartz differs from the assumptions 
that have been made so far. Quartz already undergos a 
phase transition at 847 K that has a strong influence on 
the behavior of the bands. Gervais and Piriou (|1975p found 
that in the vicinity of the phase transition the temperature 
dependence of the frequency shift for a-quartz is almost 
only related to the thermal expansion. Furthermore, the 
dampings show a linear temperature dependence up to 
700 K as would appear if only cubic anharmonicity were 
taking place, while a critical increase appears in the ap- 
proach of the phase transition. In contrast, all bands that 
survived the phase transition to /?-quartz do not show a 
strong temperature behavior any longer, even in the vicin- 
ity of the transition temperature. This behavior is simi- 
lar to that of the thermal expansion (Gervais and Piriou 
11975 ') . Another interesting point is the temperature behav- 
ior of the frequencies of the modes {E)j=2 and (^2)^=1. It 
shows a monotone increase in both TO and LO vibrations, 
whereas all other modes show a more ore less smooth de- 
crease. Both bands disappear by reaching the transition 
point. 

The possible temperature development of ec« has been 
neglected in the FPSQ calculations of Gervais and Piriou 
P974al and I1975P for corundum and quartz as well as in 
those of Thibaudeau et al. (|2006p for spinel. However, 
Cowley (|1963p has shown that Coo has a temperature de- 
pendence when taking multiphonon contributions into ac- 
count. Therefore, while fitting our data, we allowed Coo 
to vary and could indeed verify a certain temperature 
dependence for all measured materials, with the general 
trend toward an increase in Coo with increasing tempera- 
ture. This result is supported by the data of Thomas et 
al. p998p , who compared measured near IR refractive in- 
dexes of corundum with their multiphonon model. They 
found a linear increase of n with increasing temperature, 
corresponding to a quadratic increase of eoo (since the ab- 



sorption is too small to count for e in that wavelength 
range). Unfortunately, we cannot give any analytical tem- 
perature dependence of eoo until we have more data at 
higher and lower temperatures, which would increase the 
accuracy of the resulting fit. 

4. Calculation of small-particle spectra 

Using the dielectric function e (see Eq. [3]) that was calcu- 
lated during the fitting procedure, we derived the absorp- 
tion cross section for spherical particles small compared 
to the wavelength with the aid of the theory by Gustav 
Mie fMie ll90^ . If we let a be the grain radius and A the 
wavelength, then 



C-abs == — ^ — ^m 



X 



(14) 



For nonspherical particle shapes, a useful method for 
calculating small-particle spectra is the distribution of 
form factors (DFF) model fMin l20(J5|) . In the framework 
of this model, the absorption cross section is given by a 
generalization of Eq. ((T4)) in the form 



(-^abs — 



SttV 
A 



Im 



1 



P{L) 



1 + L(e-1) 3 



dL, 



(15) 



where P{L) is a distribution function of form factors L 
that are defined in the interval (0;1) and contain the in- 
formation about the particle shape. In this case, a is the 
radius of a volume-equivalent sphere. For P{L) = 5{L = 
1/3), Eq. dUl) is obtained. Mutschke et al. (2003) have 
demonstrated that certain synthetic DFFs derived by Min 
et al. (|2006p for aggregates of spherical particles and for 
Gaussian random sphere particles are well suited to re- 
producing measured 300 K absorption spectra of typical 
corundum and spinel powders with either more roundish 
or more irregular grain shapes, respectively. We calculate 
such absorption spectra at high temperatures in addition 
to those of spherical particles to demonstrate the diversity 
of spectra that can be derived from our optical constants. 
For materials with anisotropic optical behavior like 
quartz and corundum, the absorption cross section has 
to be calculated for each incoming-beam polarization par- 
allel to the crystal axes a, b, and c, and the results have 
to be summarized to an averaged absorption cross section 
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Cabs- In the case of corundum and quartz, the polariza- 
tion parallel to a and b results in identical spectra, which 
means that the absorption cross section for the ordinary 
ray {Cabso) has to be taken with double weight compared 
to the absorption cross section for the extraordinary ray 

^abs eo- 



Cabs- ^Cabseo+ ^Cabso- 



(16) 



In the case of spinel, the averaging is of course not 
necessary. It has to be noted that for all nonspherical par- 
ticle models, the separate treatment of the crystal axes in 
calculating the absorption cross sections is physically not 
exact and can lead to deviations in the predicted spectra 
from reality. This can only be circumvented by much more 
complex models such as the discrete dipole approximation. 
In the main bands of corundum and spinel, however, the 
deviations are minor, and the spectra shown in the follow- 
ing are sufficiently reliable. For details see Mutschke et al. 
((M?gi) . 

In Figs. [Mm the calculated temperature-dependent 
C'abs-spectra of corundum, spinel, and quartz are pre- 
sented. In each figure, the spectra for small spheres 
(Mie theory), aggregates of spheres and Gaussian random 
spheres are compared. For all three minerals, the differ- 
ences between the spectra of Mie theory (bottom spectra) 
and those of the DFF model (offset spectra) are the most 
prominent. In Mie theory, the strongest modes appear as 
very sharp and intense bands in the spectra of Cabs while 
already for aggregated spheres (offset -1-3), the bands be- 
come much broader toward longer wavelengths, although 
a maximum remains at the position of the Mie peak. In 
the spectra of the Gaussian random spheres (offset -1-6), 
even the maximum of the bands becomes shifted toward 
longer wavelength, while only a shoulder remains at the 
position of the Mie peak. Therefore, the FWHM of the 
bands in the spectra of Gaussian random spheres is not 
very different to the one of the aggregates of spheres. Note 
that most changes in the Cobs-spectra take place only for 
the strongest bands. The weaker bands do not show any 
strong particle shape dependence. 

5. Consequences of our measurements for the 
astronomical 13 /im band 

The small-particle spectra derived in the previous sections 
enable a better understanding of solid state bands, espe- 
cially those detected in warm and hot dust-forming astro- 
nomical environments. Indeed, the most straightforward 
application of our laboratory spectra is the onset phase 
of dust-formation in AGB stars, where mainly minerals 
belonging to the hot end of the oxygen-rich condensation 
sequence form. This is because our data refer to refractory 
minerals, and our main innovation is the systematic study 
of the high temperature spectra of these minerals. 

Interestingly, there is a subset of dust-forming stars 
on the Asymptotic Giant Branch (AGB), where refrac- 
tory dust species, such as corundum and/or spinel seem to 




15 20 

Wavelength (P-m) 



Fig. 9. Comparison of the volume normalized Cabs of 
small particles of corundum for Mie theory (multiplied by 
0.5), a DFF model for aggregated spheres (offset -1-3), and 
a DFF model for Gaussian random spheres (offset +6) at 
different temperatures. 



be over-represented compared to the silicates. (The latter 
tend to be the overwhelmingly dominating kind of dust 
for other oxygen- rich stars on the AGB.) This subset is 
characterized by optically thin shells, in most cases by a 
semiregular pulsational behavior (Sloan et al. 1996), and 
- in terms of IR spectroscopy - by the presence of the so- 
called 13 micron featureo which has been the subject of 
many papers (e.g. Glaccum 1995, Sloan et al. 1996, Posch 
et al. 1999, Fabian et al. 2001, Sloan et al. 2003, DePew 
et al. 2006), and even of a PhD thesis (DePew 2006). 

Consensus has been reached that this spectral feature 
is due to an oxide rather than to a silicate dust species. 
In fact, all three dust species discussed in this paper have 
been proposed as carrier of the 13 /xm band. Speck et al. 
(2000) propose (amorphous) Si02 as feature carrier, Posch 
et al. (1999) and Fabian et al. (2001) favor MgAl204, 
Sloan et al. (2003) and DePew et al. (2006) have strength- 
ened the case for a-Al203. It should be noted that this 
dust feature has also been detected in AGB stars outside 



^ About 90% of the semiregular variables for which IR spec- 
tra are available show the 13 /im feature according to Sloan et 



al. (1996). 
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Fig. 10. Comparison of the volume normalized Cabs of 
small particles of spinel for Mie theory (multiplied by 0.5), 
a DFF model for aggregated spheres (offset +3), and a 
DFF model for Gaussian random spheres (offset +6) at 
different temperatures. The dashed vertical line indicates 
the 13 /im band position. 



Fig. 11. Comparison of the volume normalized Cabs of 
small particles of quartz for Mie theory (multiplied by 
0.5), a DFF model for aggregated spheres (offset +3), and 
a DFF model for Gaussian random spheres (offset +6) at 
different temperatures. The features that could account 
for the astronomical 13 /zm band are actually located be- 
tween 12 and 13 /zm (hence the dashed vertical lines). 



our galaxy (e.g. Lebzelter et al. 2006) and in the spectra 
of S-type AGB stars (Smolders et al. 2012). 

As for Si02, the optical constants that we measured at 
high temperatures do not make it a more viable candidate 
carrier of the 13 /im band than the room temperature data 
did. The by far strongest band of hot Q;-Si02 is located at 
wavelength of 8.7-8.9 /im for spherical grains and also re- 
mains around 9 /im for the nonspherical grain shapes we 
examined. However, no significant narrow 8.7-9 /xm emis- 
sion feature is detected in those astronomical objects that 
show the 13 /im band. The sole (weak) resonance mode of 
hot Q;-Si02 which comes close to 13 /xm is in fact located 
around 12.5-12.7 /im, with only a weak dependence on the 
grain shape, as mentioned before and as demonstrated in 
Fig. [TT] As a result, a-Si02 cannot be considered as a 
promising candidate carrier of the 13 /im band; so we fo- 
cus on the application of our high temperature data for 
hot spinel and corundum. 

The following procedure has been applied to compare 
the small-particle spectra derived from our optical con- 
stants with astronomical spectra: i) each Cabs spectrum 
(for our measured T^ values) is multiplied with the cor- 



responding black body function B^(Td), and the result is 
normalized to unity. We refer to the resulting curves, i.e. 
to the normalized Cabs x B^(T) curves, as normalized dust 
emissivities. 

ii) Selected spectra of AGB stars showing the 13 /im bands 
are treated as outlined by Posch et al. (1999) and Fabian 
et al. (2001) to derive residual dust emissions. A polyno- 
mial is used to represent all other dust components but 
the carrier of the 13 /im band. The respective polynomial 
is then subtracted from the spectra. The result is also 
normalized to unity. A mean out of 23 such residual dust 
emissions (based on ISO spectra) was calculated. In short, 
the following comparisons confront normalized dust emis- 
sivities with mean residual dust emissions in the 11-15 /im 
range. Just for comparison, the residual emission spectrum 
of an S-type AGB star as seen with Spitzer - taken from 
Smolders et al. (2012) - is also shown. 
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5.1. Hot corundum and the 13 ^m band 

As mentioned before, crystalline a-Al203 has been con- 
sidered as a potential carrier of the 13 /im feature; e.g., 
Glaccum (1995) wrote that "hot sapphire is the most 
likely source of the 13 ^m feature found in some M and 
MS stars" (but he did not indicate on which IR data this 
assertion was based). As for the term "sapphire", we will 
avoid it in the following, since it refers to a-Al203 con- 
taining coloring impurities, e.g. small amounts of Fe^^, 
Fe^+, Ti^+ or V^+. We instead refer to pure a-Al2 03 and 
use the term corundum for it. Evidence of the formation 
of AI2O3 came from extensive research on presolar grains 
in meteorites (e.g. Nittler 1997, Clayton & Nittler 2004; 
Hoppe 2004). However, it has so far not been possible 
to set any strong constrain on the polytype of AI2O3 by 
presolar grain studies (which is a problem, since many 
polytypes of AI2O3 exist - see, e.g., Tamanai et al. 2009). 
Until very recently, it has been impossible to produce 
any satisfying fit of the 13 /im band using optical constants 
of corundum at room temperature. As shown by Posch et 
al. (1999), the following dilemma arose, based on the room 
temperature data: 

a) for small spherical particles, corundum would produce 
a 12.7 /im feature with too narrow a bandwidth; 

b) for a continuous distribution of ellipsoids, a-Al2 03 
would produce a 13.2 /im feature with an up to 7 times 
too large bandwidth. 

Takigawa et al. (2012) proposed to solving this 
dilemma by introducing ellipsoidal grains of a particular 
shape. In their condensation experiments, they found that 
the condensation rates of crystalline AI2O3 are quite dif- 
ferent for the different crystallographic axes, which leads 
to a flattening of the condensed particles along the crystal- 
lographic c-axis. They derive an axis ratio of rc/ra^O.79 
and note that for oblate spheroidal grains with a slightly 
lower ratio of rc/ra=0.70, the position of the 13 /im feature 
could be well reproduced with room temperature optical 
constants. This is because, compared to spherical grains, 
such oblate spheroids have a resonance for the E_Lc con- 
tribution at slightly longer wavelengths, while the minor 
resonance of the E||c contribution shifts to shorter wave- 
lengths. This leads to a single peak, which in the case 
of the experimentally derived axis ratio rc/rQ=0.79 would 
still be positioned a bit shortward of the 13 /im feature and 
even narrower than the spherical resonance. Allowing for 
a range of axis ratios around a somewhat changed value 
of rc/ra=0.70 would give the right feature (cf. Takigawa 
et al. 2012; see also Bohren & Huffman 1983, Sect. 5.3, 
for the applied method of calculating Cabs for ellipsoidal 
grains). We show in the Figs. |T2] and [121 how our high- 
temperature data partly confirm and partly modify these 
results. 

Figure |T2] refers only to spherical grains. The plot 
shows how the increasing temperature shifts the peak po- 
sition of the strongest resonance of spherical corundum 
grains to longer wavelengths, until a position of almost 
exactly 13 /xm is reached for temperatures close to 1000 K. 
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Fig. 12. An average profile of the astronomical 13/im 
emission band - as derived by Fabian et al. (2001) for 
a sample of 23 oxygen-rich AGB stars - compared to the 
normalized emission spectra of spherical corundum grains 
at temperatures of 300 K to 928 K. 
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Fig. 13. The case of nonspherical corundum grains. These 
grains, particularly rotational ellipsoids with an axis ra- 
tio of rc/ra=0.79, enable a good representation of the ob- 
served 13 /im band profile (the latter is the same as in the 
previous figure). In addition to our mean ISO 13 /im band 
profile, a single Spitzer-IRS 13 /im band profile (of V376 
Aur, kindly provided by K. Smolders) is also shown. 



However, the bandwidth significantly increases with tem- 
perature. At the same time, the peak height Cabs.max is 
reduced, which is not seen in the normalized spectra. It 
will be noticed that for T = 928 K, the width of the 13 /im 
corundum band is too large compared to its (average) as- 
tronomical counterpart. It should be noted, though, that 
there are individual astronomical objects with slightly 
broader 13 /im bands than others (cf. Posch et al. 1999, 
Tab. 1 and Fabian et al. 2001, Tab. 4). We also verified this 
with Spitzer-IRS spectra, kindly provided by K. Smolders. 
Broader band profiles may indicate higher (mean) dust 
temperatures. 
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Figure[T3]refers to oblate spheroidal grains with the ex- 
perimental axis ratio found by Takigawa et al. (2012), i.e. 
rc/ra=0.79. With the room temperature data, we find, in 
agreement with their result, that the calculated band is at 
a somewhat shorter wavelength and that a further flatten- 
ing of the particle shape would be necessary for a perfect 
match. For grains radiating at T=551K, however, a very 
good representation of the observed 13 /im band profile 
is achieved. For still higher temperatures, the corundum 
peak is shifted even further to the red (see the curve for 
T=738K in Fig.[T3]). This leads to a discrepancy with our 
mean (ISO-SWS-based) profile of the 13 /im band (see be- 
low). 

It should be noted that only one percent of corundum 
in the circumstellar shell of an AGB star are sufficient to 
account for the strength of the 13 /im band. This has been 
shown by DePew et al. (2006). In this study, the discrep- 
ancy between the peak position of corundum spheres and 
the observed 13 /im band had not become evident due to a 
relatively coarse wavelength grid (used to speed up the ra- 
diative transfer calculations, typically with one grid point 
every 0.5 /im around 13 /im; see their Fig. 5). 

The study by DePew et al. (2006) also reveals how it 
may be possible that we do see corundum's (or spinel's) 
~13/tm feature, while we do not see their weaker bands at 
longer wavelengths: those bands may simply be suppressed 
by temperature effects and owing to the superposition of 
other dust components (see also the final section of Sloan 
et al. 2003). 

Another argument (other than Fig. I13p for corundum 
being the carrier of the 13 /im band fact that this band is 
frequently detected together with a broad spectral feature 
peaking at '~ll-12/im, which is usually assigned to amor- 
phous AI2O3. The coexistence of crystalline corundum and 
amorphous AI2O3 in an oxygen-rich circumstellar shell 
with a low mass loss rate seems quite likely. Incomplete 
condensation of a gas of solar composition might at the 
same time result in the formation of amorphous alumina 
and in a relative lack of silicates (compared to shells with 
greater optical depths that had more time to form sili- 
cates as well, in addition to refractory oxides - see, e.g.. 
Smolders et al. 2012 for a short discussion of this idea). 

5.2. Hot spinel and the 13 fim band 

For the case of spinel, we restrict ourselves on the small 
particle spectra of spherical grains - since, as pointed 
out by DePew et al. (2006), nonspherical grains of spinel 
can hardly account for the specific properties (bandwidth, 
bandshape) of the 13 /im band. More specifically, a dis- 
tribution of ellipsoidal grains or hollow spheres leads to 
a broadening and shift to longer wavelengths of spinel's 
main MIR emission band, which is incompatible with the 
properties of the 13 /tm feature (see also Posch et al. 1999). 

In Fig. [T4j the same astronomical spectra as in the 
previous figures are shown in comparison to the high tem- 
perature emissivities of spinel grains. It becomes evident 
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Fig. 14. Comparison between the average profile of the as- 
tronomical 13/im emission band - same as in the previous 
figure - with the normalized emission spectra of spherical 
spinel grains at temperatures of 300 K to 928 K. 



from this figure that the 300 K emissivity of spinel - i.e. its 
room temperature spectrum! - comes closest to the profile 
of the 13/im band. For all the higher temperatures, the 
bandwith of spinel's strongest emission signature becomes 
much broader than for the astronomical counterpart (with 
the FWHM reaching 1.38 /xm at 928 K). The peak position 
shifts to 13.5 /im for the largest examined temperature 
928 K, which is also at odds with the peak wavelength 
of the 13/im feature. For nonspherical spinel grains, the 
13 /xm band peak position shifts to even slightly longer 
wavelengths than 13.5 /tm (see Fig. ITUl) . Table [5] summa- 
rizes essential properties of the observed 13 /tm band and 
of its potential carriers (spinel and corundum) radiating 
at different temperatures. 

Table 8. Comparison of the properties of the observed 
13 /im band profile with those of spinel and corundum at 
different temperatures, according to our new in situ high 
temperature optical constants. 



substance 


Apeak 


FWHM 


comments 


or spectrum 


[/im] 


[/tm] 




Mean ISO-SWS 


13.00 


0.45 




V376 Aur, Spitzer-IRS 


13.07 


0.60 




corundum, spher., 300 K 


12.73 


0.26 


+ shoulder 
at 13.2 /im 


corundum, spher., 551 K 


12.82 


0.35 




corundum, spher., 738 K 


12.87 


0.50 




corundum, spher., 928 K 


12.94 


0.73 




corundum, obi., 300 K 


12.92 


0.27 




corundum, obi., 551 K 


13.00 


0.37 




corundum, obi., 738 K 


13.07 


0.46 




spinel, spher., 300 K 


13.17 


0.53 




spinel, spher., 551 K 


13.35 


0.72 




spinel, spher., 738 K 


13.44 


0.97 




spinel, spher., 928 K 


13.53 


1.37 
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Fig. 15. Comparison between the average profile of the as- 
tronomical 32/im emission band with the normalized emis- 
sion spectra of spherical spinel grains at temperatures of 
300 K to 928 K. 



The peak in the mean ISO spectrum at 16.8 /im seems 
to coincide with a minor emissivity peak of spinel at room 
temperature, as noted by Posch et al. (1999) and Fabian 
et al. (2001). However, the 16.8 /xm band seen in some 
astronomical sources of the 13 /im feature has meanwhile 
been assigned to a CO2 molecular line (Sloan et al. 2003). 
Therefore, this small spectral band can no longer be used 
to support the assignment of the 13 /im band to spinel. 

There is yet another spectral feature that spinel grains, 
if present in a circumstellar shell, are supposed to produce, 
namely an emission band peaking close to 32 /xm (Fabian 
et al. 2001). Indeed, there is a narrow emission feature 
at 31.8 /Ltm in many of the sources of the 13 /im band. 
Therefore, we show an average profile of this 31.8 /im band 
(according to Posch et al. 2006) and compare it to the 
emissivity of hot spinel grains in Fig. 1151 This comparison 
shows that even for spherical grains at room temperature, 
no satisfying fit of the astronomical 31.8 /im band with 
stoichiometric spinel can be achieved (the bandwidth for 
spinel at T=300 K is too large, and the peak position is at 
too large a wavelength). For spinel at higher temperatures, 
our new optical constants lead to an even larger discrep- 
ancy between its emissivity and 31.8 /im band profile. We 
rather expect that spinel grains cooler than 300 K produce 
a 31.8 /im band with the required (narrow) bandwidth, but 
presently available data do not yet allow a decision on this 
question. 



6. Concluding remarks 

Our new optical constants show that both spinel and 
corundum grains, if radiating at a suitable narrow range 
of inermediate temperatures, can account for the average 
profile of the astronomical 13 /im band: spinel at about 
300-400 K and corundum at about 500-700 K (also de- 
pending on the grain shapes). 



At still higher grain temperatures, the resulting inte- 
grated emission, both of spinel and of corundum grains, 
becomes too broad as to fit the observed astronomical 
band. The precise temperature at which spinel's and/or 
corundum's 13 /im band becomes too broad and too 
strongly 'redshifted' compared to the observations de- 
pends on the individual astronomical object under exami- 
nation, since some observed 13 /xm bands of AGB stars are 
broader and/or peak at longer wavelengths than others. 
Therefore, the 13 /im band may even be used as a dust 
temperature indicator, at least at a somewhat more ad- 
vanced stage of astronomical MIR spectrosopy (with even 
more available MIR spectra of AGB stars and with in- 
creased signal-to- noise ratios). 

So far, we have fitted only the 13 /im band with our 
new optical constants of hot corundum and spinel. Both of 
them have additional weaker bands at longer wavelengths 
for which no counterparts in the observed astronomical 
spectra could be found so far. For spinel's 32 /im feature, 
there is an astronomical counterpart, but at a smaller peak 
wavelength and with a narrower bandwidth (only compat- 
ible with temperatures of MgAl2 04 grains below 300 K). 

The lack of (some) longer wavelength features may 
be due to a suppression of those bands by temperature 
effects, as shown by other authors (e.g. de Pew et al. 
2006). Radiative transfer calculations should be made 
again, based on our optical constants, in order to precisely 
predict the strengths of those secondary bands (at an ex- 
pected fractional abundance of corundum and/or spinel 
amounting to a few percentage points) and to check sys- 
tematically with which fitting parameters such predictions 
are indeed compatible with presently available astronom- 
ical spectra. 
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